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Abstract. The paper aims at constructing two different solutions to an elliptic system 

u-Vu + {-A)^u = XF 

defined on the two dimensional torus. It can be viewed as an elliptic regularization of the stationary Burgers 
2D system. A motivation to consider the above system comes from an examination of unusual propetries of 
the linear operator Xsin ydxW + arising from a linearization of the equation about the dominant 

part of F. We argue that the skew-symmetric part of the operator provides in some sense a smallness of 
norms of the linear operator inverse. Our analytical proof is valid for a particular force F and for A > Aq, 
m > mo sufficiently large. The main steps of the proof concern finite dimension approximation of the system 
and concentrate on analysis of features of large matrices, which resembles standard numerical analysis. Our 
analytical results are illustrated by numerical simulations. 
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1 Introduction 

Analysis of sets of solutions to elliptic systems/equations is of particular interest in Partial Differential 
Equations (PDEs) research. On one hand, the question is challenging from the viewpoint of mathematical 
techniques, on the other hand, the precise information about this set is crucial for understanding the 
dynamics of evolutionary problems standing behind the elliptic one. In general, the theory provides us with 
two answers: either there exists precisely one solution or the system admits at least one solution. 

PDEs analysis provide only few examples for quite simple problems. Here we can start from the classical 
example using the Mountain Pass Theorem for a semilinear elliptic equation |Eval0] . Nonuniqueness example 
for the stationary Navier-Stokes equations |Ga HB, important examples for geometry related to the mean 
curvature problems [BC83| . or nonuniqueness for one-dimensional viscous Burgers’ equation [BP12| . (in the 
case of evolutionary system we refer to |Dix96| . [AAlOj l. Derivation of asymptotic lower bound for the 
multiplicity of solutions for a semilinear problem can be found in [DY05] , [LM14| , and for a class of elliptic 
equations with jumping nonlinearities in |MP10] . On the other hand, the work on numerical multiplicity 
proofs for systems/higher dimensional PDEs has been an active topic of research. There exist several 
computer assisted proofs of existence of at least several solutions for particular PDEs. Let us stress that 
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contrary to our approach, all results obtained using a direct computer assistance work essentially for isolated 
parameter values or a compact set of parameter values, as all computations performed by a computer are 
finite. The results include a proof of existence of four solutions to a semilinear boundary value problem 
for particular choice of parameters |BMPn3] , an outlook for the multiplicity of solutions for some multidim 
PDEs is provided by a proof of existence of nonsymmetric solutions for a symmetric boundary value problem 
|AK12] . validated bifurcation diagrams constructed in |BLV13| . |GL11] . structure of the global attractor 
|MPMW07| . numerical existence proofs for fluid flow, and convection problems |WN09| . |HNX99j . 

The subject of the present paper is the following elliptic system, which can be viewed as an elliptic 
regularization of the stationary Burgers system |Bur48| , |Hop50| , |Col51| in 2D 

u-Vu + {-A)"^u = XF on T^. (1) 

Here u is sought as a vector function it : —)• The vector E is a particular external force, and in this 

paper we define it as 

The magnitude of the external force is controlled by the parameter A and it is assumed to be greater than 
some positive number Aq- We shall note that the system has no a-priori estimate. The issue of the existence 
of a solution to the system Q is still open for a general form of XF. To our best knowledge the case m = 1 
is even unclear. 

Let us discuss what motivated the presented research. Our numerical investigations of 0 revealed a 
solution possessing a curious structure: one of solution’s Fourier modes being of A magnitude, and the 
remainder being bounded uniformly with respect to A. We further noticed that the natural symmetry 
embedded in this equation implies the existence of a second solution, as the reflection by the symmetry 
of dominant part produces an essentially different solution. Further on, to convince ourselves that this 
structure is in fact conserved for A large values after a bifurcation, we performed a numerical bifurcation 
analysis, which showed that the graph of solution’s norm is approximately linear, and in fact there is a 
pitchfork bifurcation in the system for a particular A. 

We emphasize, that our situation is not as simple as it would be, when a symmetry embedded in the 
problem implies immediately existence of a different second solution. For sufficiently small A’s solutions 
are symmetric, and for a large A the symmetry is broken, which allows us to establish existence of two 
different solutions for large A’s. For small A’s we claim only existence of a solution, as the two solutions 
from our main result, merge into a single one. The symmetry is seen as elementary, simply enough, we can 
exchange x with y, and the first component with the second component (denoted x -H- y in the sequel). 
Apparently, for the existence of claimed solutions a stronger regularization effect, than the one provided by 
the Laplacian operator, is needed. This is why we state our main result (Theorem |1.1[) for m sufficiently 
large. Our analytical results are supported by a numerical bifurcation analysis (Section®. 

The main tool of our technique is to exploit unusual features of a linearization of the system. Let ||ic||io° 
denote the supremum norm of elements of the Fourier series w. Apparently, for the solutions to the following 
scalar problem 

AsinySa;!/; + (—A)™?/; = Asinx on T^, (3) 

we obtain 

^ ( 4 ) 

in other words, this quantity is free from A dependence (for large A), although other norms are growing 
with A. We see an interplay between the growth of the right-hand side and an increase of influence of 
the term A sin fix, which represents (in some sense) a rotation effects. In particular, it causes that the 
amplitudes of modes to be uniformly bounded. Such effect can be compared with general phenomenon of 
hypocoercivity explained in |Vil09] . We shall note, however, we do not apply the general theory for operators 
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of type A* A + B, since we want to avoid considerations in the Hilbertian spaces. We work instead in the 
framework, which is the most optimal for our analysis. The technique is elementary, yet technical, 
we perform a large matrices analysis. The features of the linear operator are indeed found for its finite 
dimensional truncation - a Galerkin approximation. The key result concerning Q is described by Theorem 


5.8| and its proof is the main part of this paper. 

Our analysis of the system 0 allows to prove the following theorem being the main result of the present 


paper. 


Theorem 1.1. Let m> 9l2 and X > Xq be sufficiently large. Then there exists at least two solutions to the 
system 0 with F given by 0 such that 

+ L^{x,y) + R^{x,y), (5) 


+ L‘^{x,y) + R^{x,y), (6) 

where L^,L^ are solutions to the linearization and they are of order in l°° norm and R^,R^ are of 

order X~°‘ with a > 0. 

The proof of Theorem o is based on a subtle analysis of the system We impose the form of 
solutions and then we construct them via approximation on finite dimensional subspaces. The natural 
symmetry x y implies that we obtain at least two different solutions, provided A is sufficiently large. 

Indeed, proved properties of the system Q are the main impact of the present paper. We are ensured 
that this type of properties will allow to study precise dynamics of systems with the transport term u -Vu. 
Here the most natural example is the Navier-Stokes equations. However, the current box of tools is not 
sufficient to attack this problem. We present here a brand new technique to study quasilinear elliptic 
systems. Hence one can look at the system Q as a toy model for which we demonstrate applicability of our 
new method. Such approaches are standard in numerical analysis and here we use them to obtain a purely 
analytical result. 

We are highly convinced that the obtained in this work explicit bounds for norms of tridiagonal differen¬ 
tial operators, which are independent of the dimension, can be applied for other problems, including bounds 
for solutions of some linear PDEs, computer assisted proofs for nonlinear PDEs, numerical analysis of dis¬ 
cretizations of certain PDEs, and slow-fast systems. There is an interest in tridiagonal operators arising in 
PDEs, see e.g. |BDL15] . Let us also note that methods based on Eourier series may be applied for systems 
in pipe-like domains. An example is |Muc03| . where analysis of the Oseen operator gave very precise space 
asympotics of solutions in front and behind an obstacle. 


and 


u^ = x( 


= X 


sm X 


Literature concerning the issue of existence to stationary Burgers equations is not rich. Most of the 
results concern only the mono-dimensional case model |BGS01[ IBP12[ IBEG"*"!!] . It motivates us to perform 
numerical analysis for various cases of the model. We observe that for the system 0 there exist a threshold 
value m > 1, such that the main result is valid for all m > m ~ the two distinct solutions can be still 
constructed. For the case m < m, especially, in case of the stationary forced 2D Burgers equations (m = 1) 
the global picture is significantly different, and for m = 1 certainly the two solutions cannot anymore be 
isolated as in the other cases. 

The paper is organized as follows. We present in Section|^the subject of this paper written in coordinates, 
in Sectionbifurcation diagrams, and a brief technical explanation. In Section]^ the relevant symmetries 
of the problem, which are crucial in our analysis. In Section the matrix form of the linearized operator, 
along with some important inverse operators bounds. In Sectiona-priori bounds for the solutions of finite 
dimensional truncations, and in Sectionj^ an existence argument for the infinite dimensional system. Finally, 
in Section some technical lemmas necessary to prove crucial inverse operators bounds from Section 
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2 Preliminaries 

We start our analysis with the preparation of our system 

u ■ Vu^ -|- = Asiny, 


u ■ + (—= Asinx. 


( 7 ) 


We fix the notation 


u = V + V, where v = X 


sm y 
0 


( 8 ) 


We focus just on construction of solution (§, the symmetry will imply existence of the second one - 
Section 4. The above relations restate the system 0 as follows 

Xsinyd^V^ + = -XcosyV"^ - V ■ VP^ 

Xsinyd^V"^ + {-A)^V^ = Asinx - 1/ • VP^. 


see 


(9) 


In order to split the solution into two parts, the first with small amplitudes and the second with higher ones, 
we introduce the linearization of Q 


A sin ydxA + (—A)™A = —A cos y B, 
X sin ydxB -|- (—A)”^B = A sin x. 


and define V as the following pair 


V = 


A 

B 


+ 


a 

b 


( 10 ) 


( 11 ) 


Then we find equations on a and b 


X sin ydxtt + (—A)”^a = —A cos yb — 
X sin ydxb + {—A)'^b = — “ + ^ 


a A 
b + B ^ 

V{b + B). 


V(a + A), 


( 12 ) 


Looking at the above problems we see that all analysis depends on the properties of the following operator 

L\{w) = Asinycla;^ -|- {—A)"^w. (13) 

The key element of the proof of Theorem |l.l| is a result concerning estimates of the inverse to Lx operator. 
The precise statement of the result we find in Section it is Theorem [HI 

Notation In bold we denote complex coefficients of vector valued functions, e.g. = (a],,a|) G C^, 
where a^, and aj, denote the first, and the second component of respectively. Let k,ki,k 2 G "L? denote 
pairs of integers. By k^, kl, k^, and kf, k"^ we denote first, and the second components respectively. 
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We write the problem 0 in Fourier’s coordinates, what is the most natural way to consider the problem 
on a torus. 


u{ 


k&I? 

F{x) = Y. Fk = {FlFl) G C^. 


(14a) 

(14b) 


fcez2 


'Y X! + (fc^)^”^) flfc — AF^ — G{a, A)-^ — 0, j — 1,2, G Z 


ki+k2=k 

ki,k2£lP 


k\+k2=k 

k-i,k2£l?‘ 


(14c) 

In order to study the problem 0 we introduce Fourier’s basis. The operator (—A)™' is diagonal in this 
basis, having eigenvalues. In order to simplify the arguments in the remainder of 

the paper we will consider (—A)™ having eigenvalues (/c^)^™ + (fc^)^™. Of course, all arguments are valid 
for the ((A:^)^ + (fe^)^)”^ case, as it bounds (fc^)^™ + (/c^)^”^ from above. 

For the particular choice of the external forcing, F in our case is given by 


-^(0,1) ~ -^^1,0) ~ b 


F,\ 


= Ff 


’(0 -1) = -^(-1 0) ~ ^k ~ ^ other cases. (Id) 

Definition 2.1. In the space of complex sequences {ak}k^id, we will say that the sequence { 0 ^} satisfies 
the reality condition iff 

Ofc = aZk, k G Z'^. (15) 

In the considered problem we impose odd periodic boundary conditions, i.e. 

u^{x,y) =-u^{-x,-y) 

u^{x,y) =u^{x + 27r,y) =u^{x,y + 2Tr) 

which on the level of the Fourier series means that we take the Fourier basis restricted to odd functions, or 
equivalently all coefficients are purely imaginary 

Re(ai) = 0, j = l,2, k£l?. (16) 

It is immediately verified that the space of coefficients satisfying (161 is invariant under the equation 0 , so 
we skip the formal calculations. Observe that (161 together with the reality condition implies automatically 
the following ’zero mass’ constraint 

j = l,2. (17) 


j = l,2, x,y £ 


ai = 0 , 


From now on we are going to consider the following finite dimensional approximation of the system (14) 


Definition 2.2. Let N > 0. We call the N-th Galerkin approximation of (14) the following system 


Y (^kfiklal^+ Y 4fik2al^ + ({k^)^"^ + {kY'^Ul-XFl = GNia,X)i = 0, j = 1,2, |fe| < A, 


ki+k2=k 

\ki\,\k2\<N 


ki+k2=k 

|fci|,|fc2|<Af 


(IP) 

Definition 2.3. We introduce the following Banach spaces for sequences {a^} equipped with the following 
norms 

||{afc}||z°° = sup \ak\, ||{afc}||zo° = sup(|fci| + |A: 2 |)^|afc|, 

fcez2 fcez2 

ii{o»)ii.. = E i“*i. ii{“'=}iii! = E (i*=ii + i*=2i)i“*i- “ 

fcez2 fcez2 

Definition 2.4. Let us define the following space 


H = 


|{ofc} G : ak = a-k, Re(afc) = 0, for 0 < |/c| < a| 


we are going to look for solutions a of (IP), such that a € H x H. In the sequel, whenever H appears, N 
will either be fixed or clear from context. 
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3 Numerical bifurcation analysis 


We analyzed the bifurcation structure of the problem (IP), we present the results on Figure]^ Starting 
from the zero solution at A = 0 we followed the branch of solutions. We detected a pitchfork bifurcation at a 


value of A depending on the m parameter in (IP I, from the point of the occurrence of pitchfork bifurcation 


we followed a stable (one of 2) and the unstable branches (it is unique). 

For a given A we were looking for a(A) such that GAr(a(A), A) = 0. We implemented the path following 
procedure in order to track a(A), for this to work the derivative is required, as the eigenvalues have 

to be controlled in order to detect bifurcation points. We implemented the path following procedure on the 
top of the software |Cyrl4| where the partial derivative is calculated by means of automatic diffrentiation 
and fast Fourier transforms. 

We present the bifurcation diagrams for two cases 

• in Figure for the choice m = 1, and varying the truncation dimension N. Our theory excludes this 
case. 


• in Figure for the choice m = 6. Our theory includes this case. 

There are some apparent differences between those two cases. In Figure and in blue we marked the 
unstable branch of index 1, and in black the stable solution(s) - this branch represents in fact two solutions 
of the same norm related with a symmetry. The symmetry is denoted by S in Section Apparently, the 
considered bifurcation is the point where S is broken. Let us relate the presented diagrams with our current 
results. In this work we prove that on the stable branch in Figure there are two distinct solutions, and 
this branch is approximately linear for any A sufficiently large. 

The diagrams were generated with the projection dimension N = 8, which corresponds to 17“^I2 degrees 
of freedom. 

To perform it we approximated the solution by a fixed number of Fouriers’ functions. On Figure we 
present a few bifurcation diagrams that were obtained using Fouriers’ approximation both of the smaller 
and larger (limited by our computational resources) approximation dimensions. To construct the diagrams, 
we started from the zero solution at A = 0, the branch of solutions (u(A)) was followed until a bifurcating 
solution was found. In case a bifurcating solution was found, both of the branches: the original, and the 
new bifurcating branch were followed. 

Observe that those diagrams significantly differ. For instance the value of A for which the numerical 
pitchfork bifurcation occurs is proportional to the approximation dimension, we mean that A was significantly 
larger, when a larger approximation dimension was used. This lead us to the conjecture that the apparent 
bifurcation is only a numerical artifact. Evidently, in the case of stationary forced 2D Burgers equations 
(m = 1) the dynamics is either not finite dimensional, or the dimension is really high. This is in contrary to 
the cases, that are included in our theory (e.g. m = 6), in which the dynamics is essentially finite dimensional 
(the bifurcation diagrams that we computed using different approximation dimensions, the same as in the 
case m = 1, does not differ much from each other). This is probably due to the not strong enough smoothing 
effect provided by the Laplacian operator compared to the higher order elliptic operators. 


4 Symmetries 

In this section we present relevant symmetries of the problem that we found, and we are going to use them 
in our arguments in the sequel. We asumme that the considered forces are also symmetric. Later on it 
will became evident that the solution we construct is symmetric with respect to S', and nonsymmetric with 
respect to S. 
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1 = 1, N=6 


1=1, N=8 




lambda 


lambda 




Figure 1: Bifurcation diagram for u ■ Vu + {—A)u = XF. Each diagram was computed with different 
approximation dimension N (given in the title). For this problem, the linear structure is not anymore 
evident, the question of the existence of two distinct solutions is left open in this case. The bifurcation point 
depends heavily on the dimension. 


6-th order Laplacian case bifurcation diagram 



Figure 2: Bifurcation diagram for u ■ Vu + (—A)®u = XF. Clearly, the stable (black) solutions curve is 
almost ||tt(A)|| = A, as there is only one mode of A magnitude. Single pitchfork bifurcation was detected at 
Ao = 8.0629. 
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Symmetry S Let S:HxH^HxH he the following symmetry (denoted x -H- y in the prequel) 

S (u^{x,y),u^{x,y)'^ = (u^{y,x),u^{y,x)^ (19) 

It is immediately verified that the solutions of the system and all its Galerkin approximations are 
invariant under this symmetry, i.e. G{a) = 0 G{Sa) = 0 as long as F is symmetric. 


Symmetry Let S^: H x H ^ H x H he the following symmetry 




( 20 ) 


It is immediately verified that the solutions of the system and all its Galerkin approximations are 
invariant under this symmetry, i.e. G{a) = 0 G{S^a) = 0 as long as F is symmetric 

Symmetry S' Let S': H x H ^ H x H he the following symmetry on the level of Fourier modes 


“)(fci,fc2) = a(_fcl,A:2)) 


{S',k^) — a(_A:i,fc 2 ) 


for k^, P e Zeven, or k^, P e Zodd, 

for k^ e Ze^en, k'^ G '^oddi or k^ G Zodd, k"^ G Zg 


( 21 ) 

It is not immediately visible that the solutions of the system 0 and all its Galerkin approximations 
are invariant under the symmetry S', i.e. G{a) = 0 G{S'a) = 0. Thus, in the next lemma we show 

that the nononlinearity is symmetric. 


Lemma 4.1. Let Nonl he the nonlinear part of (14c). Nonl satisfies 

Nonl{S'a) = S' Nonl (a). 

Proof Below, we check the symmetry of Nonfi, by the same arguments the symmetry of NonP follows. 
To verify the claim let us consider two subcases 


Case 1 k = {k^, /c^), k^, k'^ even or k^, k"^ odd. 


Nonl]^_ki^^2) - ^ + X! 

ki-\-k2=k ki-\-k2=k 




^1+^2=^ ki-hk2=k 

If we consider indices ki, k 2 such that A: = fci + A: 2 , we have that either fcj, A:| are even (odd) {k^, even 
case) or in Aj, kj one is even and the other odd (A^, odd case), j = 1,2. This implies the second equality, 
in the first term there are generated either two or none minuses from the symmetry, as both of the modes are 
from the same component, the only minus appears from the index —k^. Whereas in the second term there 
is exactly one minus generated, as the modes are from different components, this is seen clearly looking at 

(ITi. 

Case 2 k = {k^, A^), k^ even, and /c^ odd or k^ odd, and /c^ even. 


Nonlj^_ki^k2) - ®(-fci,fc2)(“^2)«(-fci,fc2) + 

ki-\-k2=k ki-\-k2=k 


-kl,kf)^2al-kl,q) 


+ E -(‘5'“)(fcbA:2)(-fc2)(^'«)(fcl,A:2)+ E (‘^'“)(fc;,fc2)/ci(5'a)[^i^^2) = iVon/f;-. 

ki-\-k2=k ki-\-k2=k 
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If we consider indices fci, k 2 such that k = ki + k 2 , we have that both of the indices in one of the pairs 
fej, A:| are even (odd), and in the second pair indices are of different oddness (one even, and the other odd). 
This implies the last equality, in the first term there is one minus generated from the symmetry, as both of 
the modes are from the same component, the second minus appears from the index —k^- Whereas in the 
second term, as the modes are from different components, there are either two minuses generated {k\ even, 
ki odd or vice-versa, and fel even (odd)), or none minuses (/c(, fef even (odd), and k^ even, odd or 
vice-versa). 

□ 


5 Linearized finite dimensional system (truncated operator) 

We present the form of the linear operator 

Lx{w) = \s\nydxW + (22) 

in Fouriers’ coordinates introduced previously. Here we argue how to reduce the problem of deriving dimen¬ 
sion independent bounds for the inverse of (22| to the problem of bounding particular matrix norms. The 
operator (—A)™" in coordinates is diagonal 

(-A)"* = diag(...,/2"* +A:^™,...). 

To show the action of the \s\nydxW component, we introduce the following subspaces 

Definition 5.1. Let I G {0,1,..., iV}. We denote 

H D = {a: a ^ span {sin(/x), sin(±/x -|- ky) for A: = 1,2,..., N}} . 

It is easy to see that subspaces are invariant for the operator L\{w), i.e. L^H^ C HK And moreover, 
H can be decomposed 

H © • • • © and H = H, 


Hn)- 


(23) 


Definition 5.2. Let us denote the projeetion of Lx onto by 

L\ :=PhiLx: ^ hK 

Define the projection of Lx onto the following space 
Definition 5.3. Let R C H denote the following spaee 

R:= 

The projeetion of Lx onto R will be denoted by 

Lx := PrLx : R ^ R. 

To understand the action of Lx on a vector from we take 3 a = oq sin lx + sin Plx + jy 


Xsinyox + Xie 


^ ^ P c.c. 


(24) 


k=l 


k=l 


, / N-1 iV+1 

^ \ k=l k=l 




+ C.C.] , 
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where c.c. denotes complex conjugates. Therefore 

= -Mat - Ma^ + (25a) 

(Lia)^ . = Ma+_, - Ma+^, + (P™ + f^aj, (25b) 

(L^a) . = -Ma-_^ + MaJ^^ + + f^)aj , for j > 0. (25c) 

Definition 5.4. Let H’ he the following spaee 

H' = {ae H X H: a satisfies S'a = a, and S^a = a} . 

From now on, we look for solutions a G H', i.e. symmetric with respect to the symmetries S', and S^. 
Every time the operator Lx appears, it is assumed to be acting on a component of a G H'. Now, let us 
analyze the action of Lx on each of the components of a G H'. Note that for such solutions we may take 
the projection of Lx onto the subspace 


span {sin(/x), sin(/x + ky) for A: = 1,2,..., N} , 


as the modes a~ appearing in ( [25| ) are related with the modes aj through 5^ (we have S'a = a). We 
deduce from ( |25[ ) that Lx acting on the projection of each of the components oi a £ H' onto the subspace 
span{sin(Zx),sin(/x + ky)} has the following tridiagonal matrix form, which we are going to exploit in the 
sequel 


— 

— 


l2m 

IX 

0 

0 

0 


-IX 0 
2™ + 1 -IX 


IX /2m + (iV _ l)2m 
0 IX 


0 

0 

0 

-IX 

j^2m _j_ 


(26) 


The full linear operator Lx acting on each component oi a £ H' has the following block diagonal form 


Lx = 


^x 


0 

0 


0 

Li 


0 L 


N-l 

X 


tN 

^X 


(27) 


Bounds for matrices inverse to L^, Lx 

In this part we provide results on bounds of the particular norms of inverse tridiagonal matrices. Some 
technical lemmas used to prove the presented bounds are provided in section 

Lemma 5.5. Let N >0, m > 1, I = 1,..., N. The following uniform bound holds 


{Li 


^2m 

< fork,j = l,. 


,N. 


Lemma 5.6. Let N >0, m > 1, I = 1,..., N. There exist C{m) > 0 (independent of A and N), such that 
for A > 1 the following bounds hold 

E \{Li~\\<C{m){lX)-^+^/^^, E \{Li~%\<C{m){lX)-^+^/^'^, 

j=l,...,N j=l,...,N 


for all i = 1,..., N. 
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composed of bounds for the 
Where first two are standard 

norms, and the third (which we call the gradient norm) is defined 


Definition 5.7. Let A be a block-diagonal matrix 




' Ai 0 

0 

0 

A = 

0 0 

An-1 

0 


0 ••• 

0 

An _ 


where Ai,^ Ajq are N dimensional square matriees. 

We eall the gradient norm the following matrix norm 


In the next theorem we present the main result of this section, whi ch i 
following norms \\L^~ ||;i^p, \\L^~ \\L^~ see Definition 


2.3 


= . max ||A(j)||n 


max ||^z(?')||;i 
1=1,...,n" 

j=l,...,N 


max 
1=1,...,N 
j=l,...,N 


N 

'^{l + k) 


k=l 


(A) 


k,j 


(28) 


where A(j) denotes the j-th column of A. 

Theorem 5.8. Let I = 1,..., N. Let be the matrix given by (26), L^ be the truneated matrix PrL\. 
The following estimates hold for the matriees ^ (diagonal submatriees of L^ ^). 


- 1 , 


\L\ < Ci(m)((A)-i+V2m^ 




< 2‘^^{IX) 


-1 


IT 


I -i| 


,1^,! < C’2(m)(fA)-‘+'''’". 


The following estimates hold for the matrix L^ ^ 

II r^~^ll ^ r)2m\—1 

llTf'llzi^zi <Cl2(m)A-i+V-. 
We present a proof of this theorem in Section 


6 Fixed point argument 


Having the estimate for the operator (T^) ^ we are prepared to prove the main result of the paper. Con¬ 
sidered solutions to are assumed to be finite dimensional. It allows to use the results for the matrices 
presented in Section Let us define two projections of the space H 

Pr projection onto © • • • © (the rotation like part) , (29) 

Pd = L — Pr (the diagonal part) , (30) 


where Pr is the projection onto the subspace free of y dependence. We proceed as follows. First, we construct 
an a-priori estimate for the solution of (10). Let us display basic features of {A,B) € H' - solutions to (10), 
which follows directly from the bounds presented in Theorem 5.8 Observe as PrA, PrB = 0 it is enough 
to use the bound for (L^)“^, and we obtain 


^llo < < aV™, 

B\\ii < AV- P||,i < A3/2-, 


(31) 
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Consequently 


\\{A,Bf-VAy < ||(Ai?f||ziP||,j < A1/-A3/2-< A5/2-, (32) 

\\{A,Bf-VBy < ||(Ai?)^lbip|Li < aV-aV- < a2/-, (33) 

\\{A,Bf ■ V^lboc < \\{A,Bf\\io.\\A^ < aV2-a3/2- < A^/™. (34) 

In the estimations above, and generally in the estimations performed in this section we use often Young’s 

inequality for products, i.e. ||/*5lli < ll/l|i||ffl|i, and ||/*5lloo < ||/||oo||5lli- 
Now we denote V solution to the system (|^ in the following way 



where {A, By is the solution of the linearized system (10). 

First, we present the crucial property, which allowed us to establish the fixed point argument presented 
in this section. We look for solutions such that (a, 6)^ G H'^ and it follows from the symmetry S^{ax, by)'^ = 
{by, ttx)'^ that we have 


abx + bby = abx + bax = {ab)^ , and in consequence Fob = 0. 


(36) 


Standing assumptions. At the formal level of the a priori estimate we assume that solutions to (12) 
fulfill 


{a,bf <Ai-V-, 

IIFroIIp < 1, 

m > 9/2. 


Recall (12) 


A sin ydxtt + (—A)”*a = —A cos yb — 


A sin ydxb + (—A)™'6 = — 


d A 
b + B 


a A \ 
b + B )• 

V{b + B). 


V(a + A), 


(37a) 

(37b) 

(37c) 


(38) 


In order to find the bound we apply formally the estimates for L^, assuming that the solutions are finite 
dimensional. To get the bound we treat the rhs of (38), and we have 


Bound for ||6||,i = ||Pr6||;i. 


Il^ll«l ^ ll-^A ILw/l (ll(^) + II(A, B)^ • V6||p^ 

< A-'+i/” (a"/” + ||<i||p||6||,. + I|6||p||6||,. + llallpA'/” + IlftllpA'/” + A>/”'||(,||,.) 

< A-l+3/m^||^||^iA-l+2/m, ( 39 ^ 


where the last inequality is after cleaning the absorbed terms, which is due to the assumptions (37a), and 
(37c). We will also need the estimate for ||6||p, derived analogously as above 




a||pA-^+=^/2m 


(40) 


Let us define 

WPoaWi^^ = sup A;2'"|a(o,fc)|. 
|A:|<Af 


(41) 
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Bound for ||PDa||j°° . For this case the operator PdL\ is diagonal, therefore we bound the particular 
norm ||-Pz)a||«“ , it is trivially bounded by the norm of the right hand side. Moreover, observe that 
norm bounds i.e. we have ||PDa||;i < ||PDa||;|^ for m > 3/2, remembering that the dimension is two. 

^ ' '^^11*°° + \\Pd{clO,x)\\i°° + \\PDibay)\\i^ + ||z°° + II-Pd(^^j/) l|z°° 

+ ||PD({ABg-Va)||,« 

< A^/"* + ||PRa||,.||P„a||,. + ||6||,.||P„o||,. + ||P„a||,.A3/2"* + ||6 ||,.a!>/2'" + A‘/’"||Pro||,.. (42) 


We removed all terms, which do not generate Pd, i.e. any product of terms in which one is Pd, the other 
Pr, and when the bound (42) is used (potentially the worst term \\PDaPDax\\ is not present as Pd^x = 0) 
we get 


IIPd»IIi?., W™ + IIPRollii + A-i+!'/’"||Ph<.||,, + ||<.||,.||Pro||,.A-m/"' 

+ ||Pro||pA3/2™ + a- 1+9/2"^ + ||a||pA-i+^/2m ^ x^/m\\Pj^a\\ii 

<A2/m_ 


To get last inequality we used the assumption (37b I and (37c), the term A^/"* is clearly of the highest 
order from the terms that are left. Here we use that m > 9/2. 


Bound for ||PRa||;i. Observe that we have ||PDa||;i < UPdoII/^ ^ A^/™' 

II-Pr«IIzi i+i/m ||a||6||p + II(H, B)^ ■ VH||;i + ||a||p ||o||;i + ||&||/i ||a||ii + ||o||;i ||H||;i 

+ ll^llzi ll^llzi + ll(^!ll«llzi) 

<A-‘+>/'” (a|| 6||,, + A=/"” + IIPRall?, + ||Pr<i||,, A"/” + ||6||,, ||PRa||,| + 

y2/myZ/2m ^ || p^o||;l ^ ^l/my2/m\^ 

<A-‘+>/”> (a|| 6||,. + ||P„a||2. + \\PRa\\,,\^l’^ + + ||6||,. ||P«a||,. + A^^”) . 


Observe that after the second inequality the term ||ft||p||a||p is not present as Pdo-x = 0, clearly the highest 
order term is ||Pdo||p ||H||;i = 

Now we use the bound (|42|), and remove some of the terms that were absorbed by using the assumptions 


(37a), (37b), and (37c), observe in the inequality above the bad looking term A||6||p, we estimate it using 

(dof 

A||6||p < A5/2- + a 3/2™ (||P^a||p + llPfiallp), 


WPnaWil < A-i+V- (a 5/2- + A2/-A3/2™ + ||P«a||| + ||PKa||;iA2/-+ 

After using the assumption (37b) all terms with ||P/ja||;i are being absorbed, and clearly the highest order 
term in the parenthesis is so finally we end up with 


I|Pr»II, 1 < A-'+i’/"”. 


(43) 


Observe that Pru is mapped into itself by the operator due to the assumption (37c), namely m > 9/2. 
Going back to (39) we get that 

||6||,. < A-‘+*/'". (44) 

Summing up the considerations from this part we obtain the following result 
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Lemma 6.1. Let a,b be a small solution to problem (381, then it obeys the following dimension independent 
a-priori estimate 


< c-A-1+4/-, \\PRa\\n < CA-1+9/2. 


WPnaWii < 


(45) 


7 Proof of main theorem 

Using the so far presented results, we may now proceed to proving our main result - Theorem |1.1[ Here we 


want to construct the solutions, using the system (38) and the a-priori estimates (Lemma 

We start with the construction of the sequence of solution’s approximations. We define the solution 
(on+i, bn+i) as the solution to the following problem 


Xsinyd^On+i + (-A)™-a„+i = -Xcosybn+i - 

Oin T \ 


Xsm.ydxbn+1 + {-iX)'^bn+i = - 


bn T Pr 


( a-n ~\~ An \ 

+ Bn j 

V {bn -|- Bn)- 


V {On + An), 


(46) 


We take (ao,6o) = (0,0) and define An, Bn as the projection onto the spaces H(jv+ 2 )) where N determines 
the number of active modes. If {an,bn) G H(m, then {an+i,bn+i) G H{ 2 N+ 2 )- 
Repeating the estimates for the system (381 we find that if 


ll^nlLl < CX-^+^/^, WPRanWn < CA-l+9/2m^ 


WPoanWii ^ 


< CA^/"*. 


then 


||6„+i||;i < CA-i+^/™, ||PKa„+i||;i < CA-^+9/2™, WPoan+iWii < CA^/^ 


(47) 


(48) 


with the same constants C, provided A sufficiently large. In addition the form of the system (38) preserves 
the symmetry S^, too. 

We shall underline that for a fixed n we are allowed to apply results for the finite dimensional approxi¬ 


mation of Lx- We emphasize that all constants in Theorem 5.8 are N independent. 

We want to prove that {on, bn} is a Cauchy sequence. We consider the following system 

A sin yd^{an+i - a„) -h (-A)™(an+i - On) = -X cosy (bn+i - bn) 


O-n T An 
bn T Bn 


V(Un + An) -|- 


O-n—l T An—1 

bn-1 + Bn-1 


Xsmyd^{bn+i - bn) + (-A)™(6„+i - bn) = 


0-n T An 
bn T B-n 


V {bn + Bn) -|- 


0-n-l + An-l 


XX{On—l T l)) 


V {bn-1 + Bn-l)- 


(49) 


^ ^ ^ bn—l T Bn—1 

Taking a large n we want to prove that 

||RR(®n-|-l 0,n),bn+l T A ^ ||H£)(cijj-|_i On 

2 (ll^it(«n - 

where e„ —>■ 0 as n —oo, the quantity Cn is related by norms of terms like {An — An-i) and {Bn — Bn-i)- 




1 < 
1 


l),bn — bn-lWij + X '^^"^\\BD{o-n — CLn-l)\\lj^ + £n, (50) 


In order to justify (|50|) we point out few estimates which provides the inequality. Here we use the same 

). We have to estimate the following terms. 


tools as for proving Lemma 6.1 Hence we estimate the rhs of 


For \\bn+i — bnWp we have 


O-n T An 
bn T Bn 


XX{Bn — Bn-l)\\ll < Cn- 


(51) 
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( 52 ) 


For n sufficiently large it it clear that \\Bn — Bn-i\\ii —)■ 0 as n —)■ oo. Next, 

and 



ttn—1 ~\~ An An—l 
^n—1 H“ Bn Bn—I 


* V(6ri-1 + -Bn-l)!!;! 


< + A'/”^)||(a„ - a„_i),6„ - 6„_i||;i 

< e„ + (A-1+"/'" + AV-)||P^(a„ - a„_i), - 6„_i||;i 

^ (^-i+4/m ^ A^/'")A2/"*A-2/™||PB(an 


«n-l)||/i. (53) 


Hence 


\\bn+l - bnWil - - an-l),bn - 6„_l||;l 

+ (a- 2+7/™ + A-'+^/”^)A-2/-||P^(a„ - a„_i)||,i + e„. (54) 

For \\PD{an+i - an)\\ii we have 

lift (( 11 b” ) ■ ~ - “”-■>11'!' (“) 

The remaining terms here are simpler. So 

WPoian+i - On) 11^1 < A^/^™'||Pr(o„ - On-i),6n “ bn-iWii + better terms. (57) 

For ||PR(an+i - On) 11^1 we have 


ll^« {(ll + il)- - «-i) ) 11^^ ^ A2/-||P«(an-On-i)|Li+(A-i+"/-+AV-)A2/™A-2/-||P^(an-On-i)|Li 


and 


\\Pr[{ 1 11^1 + (59) 


The last term is ||Acosy(6n+i — bn)\\ii, using the estimates for ||L^||p^;i from Theorem 


5.8 


we find 


||Acosy(6n+i - bn)\\ii < X^/^n\RHSm 


^l/2m (A2/™||PH(On - an-l),6n " 6n-l|Ll 

+A3/-A-2/-||PB(an - On-l)||zi + Cn) . (60) 


(58) 
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Hence 

\\PR{an+i - an)\\ii (61) 

(A5/2-||p^(a„ - an-i),bn - bn-iWii + ||(a„ - an-i)\\ii + 

(62) 

<A-^+"/2™||P«(a„ - an-i),bn - bn-i\\ii + A'^/™||(a^ - a„_i)|bi + e„. (63) 

Summing up (54|,(57) and we conclude 

||-PR(^n+l 0 ) 77 ,), ^72+1 ^n||/i ^ ^ ||-PD(o^n l) ||/i ^ 

A-'+"/"™||PMan-a„_l),6„-6„_l||n (64) 


For m> 9l2 and large A we got (50). 

The condition (|50|) implies that the sequence {PRan,hn,PDO.n} has a limit in the space l\. It means 


that there exists a solution to problem (38) obeying estimates from Lemma 6.1 In other words we have 
constructed the solution (|^. We shall underline that the limit in l\ implies that the derivative is uniformly 
bounded, thus the nonlinear term is described pointwisely. A bootstrap method implies that the solutions 
constructed in the above way are indeed smooth. Existence of the solution (^ follows from the symmetry 
X -H- y. Largeness of A implies there are two different solutions. Theorem [LT is proved. 


8 Technical lemmas and proof of Theorem 5.8 


Notation Let A > 0 be an even number, m > 1, 1 > 0, A G 
Let us denote 


d2x2 


3T\a,b) = 


a —IX 
IX b 


We denote a tridiagonal matrix with elements {ajjjLi on the diagonal, —IX over diagonal, and IX under 
diagonal by 


iNxN 


3 T*(ai,..., Oat) = 


ai —IX 0 0 0 
IX cl2 —IX 0 0 


0 IX ojv 

d[ = 4 = /2m ^ ^ fm + _ l)2m 


Let the increasing sequence be given by 


\N 

h= 

;2m 1 / 


We denote the tridiagonal matrix with the increasing sequence on the diagonal by 

Lemma 8.1. Let I > 1. Let the sequenees { 0 ^}, {b^j} be given by the following reeursive formulas 


“0 ~ 0) 


4 = 




N 


4-i4 + ^"a2’ 
I d^N-2 + 


6[) — 0, 


b[ = 


d[ 


+ /2A2’ 


®2 ~ 


J 


4_34_2 + ai4-3^^A2 + /2A2 ^ 
d(v-2i+2 + 4-1^^A2 


al ul /2\2_|_/2\2' 

“Ar-2j+l“iV-2j+2 ^ A ^ L A 


b\ = 


6 ■ = 


4 + h{PX^ 


44 + 6i4l2A2+/2A2’ 

4j-i + b]_^l^X^ 

4,4-l + ^tWU'A2 + /2A2’ 
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N=50, m=2 


a_k + 

-^+11 b.k X - 

. + -t- ^ 'y=b.001 - 


X X y 
X + ^ ^ X , 




0.0004 - + 


0 - X X • 



N=50, m=4 



N=500, m=4 



Figure 3: Graph showing numerically calculated redfisive series {oj}, {bj} for A = 1000, and for m = 2 
(top figures), m = 4 (bottom figures). Apparent upper bound for the series is also shown (1/A). 












for j > 1. Then the following bounds hold for all j > 0 

22m 2 

0<b‘<-. 

Proof The Oj, Uj > 0 part of the bound is trivial. 

Now we prove 2^™'(ZA)~^ > We proceed by induction, first we prove that 2^™'(ZA)“^ > a\ holds. 
Observe that for all > 1, and / > 1 we have 

l2m + (jv - 1)2™ 


^ /2m 22™-1 [{N - 2)2™ + 1] ^ 

^ /.O- /HT ^ \ n - \ 0 .0.0 ^ 


(/2™ + (AT - l)2m)(/2m + (/V _ 2)2™) + ^pm + (/V _ 2)2™)2 + p\2 


~)2m |^^2m _j_ 


■)2m. 


i2m 


■)2m 


< 


(P™ + (IV - 2)2™)2 +/2A2 'a2 + /2A2 ^ ' /A 


(65) 


where a = + {N — 2)^"^, we used the estimate due to convexity {N — 1)2™ < 2 ^™ 


2m 


< 


■)2m—1 


{{N — 2)^”^ + 1 ), the last inequality follows from + 1“^)? > alX. 


Assuming 2^™(/A) ^ > aj-i we verify that 2^™'(^A) ^ > aj holds. 

d} 

First, observe that /(a) := ji -- +Px 2 is a strictly increasing function for all a > 0 


(denominator is positive), as 

/(a) = 

so we have 

4+2 + aj-il^X‘^ 


^N-2j + 2^N-2j + l 


-^N-2j + l 

l^X^ 


{^N-2j+2(^N-2j+l + ^^‘^^‘^'^N-2j+l + 

/2™ + (A:+1)2™ + 22™/A 


> 0 , 


( 66 ) 


')2m 


< 


4+24+1 + «i-l4+l^^i^^ + ~ + {k + l)2™)(/2m ^2m) 22™/A(P™ + P™) + PX^ 


^k+1 


^ IT- 


where k = N — 2j. The last inequality reduces to 

pm 22 ™-l(/c 2 m + 1) + 2‘^^IX 


')2m 


< 


(/2™ + A:2™)2 + 22™/A(/2™ + A:2™) + ^2^2 - /A ’ 

after grouping the terms in this inequality it is easy to see that it is satisfied for all / > 1, and A: > 0. 

Obviously {IX)~^ > bi holds. Analogically as above, assuming (/A)“^ > bj we verify that {IX)~^ > bj+i 
holds (it can also be verified that ffb) is strictly positive, and it is enough to verify the inequality setting 

b, = (/A)-i). 

d2i+l + bj.iPX^ ^ /2m ^ (2j)2™ + IX 


< 


1 

< — . 


d2j+2d2j+l + bj-id2j+2PX^ + PX^ - (/2™ + (2j + 1)2™)(/2™ + (2j)2™) + /A(/2™ + (2j + 1)2™) + PX2 - /A’ 

( 68 ) 

The last inequality holds due to following inequality, which is clearly satisfied 


a T IX ^ 1 
z2 + alX + /2A2 - /A’ 


where a = + (2j) 


:\2m 


□ 


Lemma 8.2. Let I > 1. All elements in 2 x 2 diagonal bloeks of L\ ^ are estimated uniformly. Preeisely, 
the following inequalities hold for j = 0,, N/2 — 1 


^ 22^771 

\iL{ )2j+l,2i+l| < 


\{L{ )2i+l,2i+2| < 

\{P\ )2i+2,2j+l| < \{L{ )2j+2,2j+2\ < 
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Proof Here we assume that I is fixed, and we drop the superscript in the notation of T\ 0 }^., 6^, and dj,, 
we use simply T, a^,, b^, and dk respectively. Moreover we denote A = IX. First, all matrices considered are 
invertible, which is obvious by calculating the determinant of tridiagonal matrices. 


We are going to use the following convention for block decomposition of 


i -1 


tI -1 _ 
— 


ji - 
^x — 


W(di,.. 

■, dAr)ii 

kF(di,. 



r(di,..., dN- 2 ) 

-A 

W(di,.. 

•, dAr)2i 

kF(di,. 

• , dAr)22 



T{dN-i,dN) 

■ V{dn. 

• • , dAr)ii 

V{dn. 



' T{dnd 2 ) 

-A'^ 

V{dn. 

to 

V{dn. 

• , dAr)22 


A' T(d3, 

1 


-1 


-1 


(69a) 

(69b) 


where 


■ 0 

0 ■ 


■ 0 

A ■ 

0 

0 

, ^' = 

0 

0 

_ A 

0 _ 


_ 0 

0 _ 


A = 


We will call Wj^k, yj,k the inverse blocks. The explicit formulas for the inverse blocks are obtained from 
the following system of equations (simplifying the notation by dropping the brackets with parameters, i.e. 
W = W{di,, dTv), Ti = r(di,..., dN- 2 ), Tjj = T{dN-i,dN) etc.) 


TjWn + AW 21 = I, 

-A^Wn + T11W21 = 0 , 


T1W12 + AW22 = 0 , 

-A^Wi2 + r//W22 = /. 


(70a) 

(70b) 


When the equations for diagonal blocks are decoupled we obtain 


Wn = 
W 22 = 


T(di,... ,djv-2) + AT(djv-i,dM) ^A?' 


n-l 


where 


T{dN-i,dN)+A T{di,... ,dN-2) A 

. 0 

, 0 


n-l 


AT{dN-i,dN) A = 


0 ... 0 0 

[0 ... 0 A^r(dAr_i, d7v)i/ 


A^ r(di,...,djv- 2 )“M = 




0 


A^T{di, ..., dj\i-2) 

0 0 


Now, we state the crucial observation - the inverse diagonal blocks Wn and H 22 are inverses of tridiagonal 
matrices, i.e. 


VFii = r((ii, ..., dAT-s, dAr_2 + A^ai) 
W22 = T{dN-i + X^T{di, ..., dN-2)~MM-’ 


(71a) 

(71b) 


where ai = T{dN-i,dM)n = 
symmetric calculations, i.e. 


An 


dN-\dN+X^ 


The same holds for the diagonal inverse blocks Vn and V 22 by 


Fii = 

F22 = 


r(di,d2) + H'^r(d3,...,d^)-M'] ' = r(di, d2 + A2T(d3,...,d^)r/) 


r(d3,...,d,v)+^'T(di,d2)-M'^ ^ =r(d3 + A26i, d4,...,djv) 


-1 




(72a) 

(72b) 
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where = T{dud 2 )^i = 

Observe that the decoupling of the diagonal blocks described above can be iterated, and the matrix 
W{di,..., d 7 v)ii = T{di ,..., dN- 3 , dN -2 + is further decomposed 


r((ii,..., d7v-3, c?iV-2 + A^ai) ^ 


W{di,..., dN-3, dN-2 + A^ai)ii W{di,..., dN-3, dN-2 + A^ai)i2 
W{di,..., dN-3, dN-2 + A^ai)2i W{di,..., dN-3, dN-2 + A^ai)22 

T{di,... ,dN-4) —A 

T{dN-3,dN-2 + X^ai) 


thus (simplifying W := W{di ,..., dN- 3 , dN -2 + A^oi)), we write the formula for the inverse diagonal block 


VFii 


T{di ,..., dN-i) + AT{dN-3, dN-2 + A^oi) ^A^ 


r(di,..., (ijv_5, djv_4 + A^a2) 


where 


a2 = T{dN-3,dN-2 += 


dN-2 + UlA^ 


dN-3dN-2 + aidN-3^^ + A^ 

From repeating j times the procedure of taking the upper-left inverse diagonal block and decompose it 
further like in (69a I, we obtain the explicit formula for the N — 2j dimensional upper-left inverse diagonal 

block of L\ ^ (below we denote Wn = W (^di, • • • , (i 7 v- 2 j-i-ij dAr_ 2 j +2 + ) 


= 


T{di ,..., dN-2j) + AT{dN-2j+i-,dN-2j+2 + A^Oj-i) 


-1 


= T{di, . . . , dN-2j-l, dN-2j + A Qj) 


-1 


(73) 

(74) 


where 


ij — TiydN—2j+^i dn—2j+2 4“ A — 


dN-2j+2 + Ojr-lA^ 


dAr-2j-|-ld7V-2i+2 + aj-ldAr_2j+lA^ -I- 

Performing iteratively j times the symmetric procedure to the one described above (performing decompo¬ 
sition like in ( 69b[ )), we obtain the explicit formula for N — 2j dimensional lower right diagonal inverse 
block 


4^22 


T{d2j, ■ ■ ■ ,dN) + A'T{d2j-2 + A^6j_i, d2j-i) ^A’'^ 
T{d2j + ^^bj, d2j+i, ..., dN)~^, 


(75) 

(76) 


where 


bj = T{d2j-2 + A^6j_l, d2j-l)22^ = 


d2-i- 


2j-2 + 


d2j-2d2j-i + bj-id2j-iX‘^ + A^ 


Note that the recursive series {aj}, {bj} are generated from the procedures described above. Using above 

results, we may now derive an explicit formulas for the 2x2 diagonal blocks of Let us present an 

example how it is done 
z-i 


Observe that from (731 we have that the N — 2j dimensional upper-left block of 
equals to T{di ,..., dN- 2 j-i, dN- 2 j + A^aj)“^ 


Then, we apply to T{di ,..., dN- 2 j-i,dN- 2 j + A^Oj) 


-1 


L 

j times the procedure of taking the upper left diagonal block, and decomposing like in (69bI, and we get 
that the j-th (counting from the bottom) 2x2 diagonal block of equals to 


T{dN-2j-i + A bk,dN-2j + A aj) = 


dN-2j-i + A^6fc —A 

A dN-2j + Attj 


n -1 


where k = {N — 2j — 2)/2. 


20 
























Let us denote D = dj^- 2 j-idN- 2 j + we have 

T{dN-2j-i + A^6fc, dN-2j + ^ —-, 

T{dN-2j-i + A^6fc, dN-2j + A^«i)22^ ~ ^ D —~— 

T{dN-2j-i + A^^fc, div-2j + ^‘^aj)i2 _=*=A 

T{dN-2j-l + A^5fc,dAr_2j + ^‘^0'j)21 dd 


We have that for f{aj, bk) = 


dN 23 +^ 0-3 _ N partial derivatives equal to 

df{aj,bk) -X^N^ 


df{aj,bk) _ A^ „ 

dai ’ 


dbk 


D2 


< 0 . 


Hence, to bound T(djv- 2 j-i+A^^fc 5 dN-2j+X‘^0'j)u upper end of the bound for aj from Lemma 

i.e. we se' 

bounding 


•)2m\ — l ^ 

diV-2j+2^’"A 


8.1 


i.e. we set aj = 2‘^^X and we use the lower end of the bound for 6^, i.e. we set b^ = 0. We are left with 


A2+djv_2j-irfiv-2:7+'ijv-2j-i2^™A 

To bound T(dN- 2 j-i + X^bk, dN- 2 j + X^aj) 


are left with bounding 


dN-2j-l+^ 


A^+djv _2j-lrfjV-22+djV-2j A 


which was already showed in (67) to be bounded by 2^”^A 

analogously as above, we set aj = 0, and bk = A“^, and we 
which was already showed in (68) to be bounded by A“^. To 


-1 
22 > 


bound the remaining two elements, i.e. T{dN- 2 j-i+X^bk, dN- 2 j+X^o,j)i 2 , T{dN- 2 j-i+X^bk, dN- 2 j+X^aj) 2 i 
we set Oj = 0, = 0, and we obtain the claimed bounds immediately. 


□ 


Lemma 8.3. Let I > 0. The following uniform bound hold 

|(Lr|)2,+l,2.+l| < ^ |(Lr|)2,+l,2fc+2| < ^ 

\{L{ )2i+2,2fc+l| < TX K-^A )2i+2,2fc+2| < /y 


fork,j = 0,...,[{N-l)/2]. 


(77a) 


Proof Here we assume that I is fixed, and we drop the superscript in the notation of T\ a^, 6(,, and 
we use simply T, a^, bk, and dk respectively. Moreover we denote A = IX. 

First, for the sake of presentation, let us prove that the claimed bounds are true for the 4x4 upper left 
corner submatrix of ) i-®- T{di, d 2 , ds, d^ + A^ajv/ 2 - 2 )~^) the general result will follow 


W (di, d2, ds, d4 + A^a^/2-2)ii tT (di, d2, ds, d4 + A^a7\r/2-2)i2 
W (di, d2, da, d4 + A^a7v/2-2)2i W (di, d2, da, d^ + A^ajv/2-2)22 


r(di,d2) -A 

AJ' T{d^, di + A^a7v/2-2) 


From the equations for inverse blocks (70) it follows that the block beyond diagonal satisfies (denoting 
W := VF(di,d2,d3,d4 + A^aAr/ 2 - 2 )) 


W 21 = r(da, di + X^a^/ 2 - 2 )~^A^Wii. 


8.1 


From Lemma 
A^aAr/ 2 - 2 )“^ are t 


follows that 2^'^X ^ > ajv/ 2-2 > 0) hence the bounds for all elements of T{d^,di + 
re same as those derived in Lemma 


1.2 


Observe that 


W 21 = r(da,d4 + X\i,/ 2 - 2 )-^A^Wii = 

T{ds,di + A^a7v/2-2)i/A(lTii)2i T{d^,di + A^ayr/2-2)i/A(VFii)22 
r(da, di + X^aj^i2-2)2i A(kFii)2i T(da, di + A^a7v/2-2)2/A(Wii)22 


21 

























as we have from Lemma 
:2„ . ^-l 


8.2 


the bounds |A(Tyii) 2 i|, |A(VLii )22 
\nd3,d^ + X^a^/2-2)u\<2^^X-^ 


satisfy the following bounds 


\T{d3,d4 + X a]\f/2-2)2i 


I < 1, and 

< A“^ . Elements from the block IL 21 clearly 


|(VL2i)ii|, |(VF2i)i2| < 2^™’A |(1L2i)2i|, 1 ( 1 ^ 21 ) 22 ! < A 


The block IT 12 satisfies symmetric bounds by a symmetric argument. Observe that the bounds for the 
diagonal blocks Wn, 1^22 were derived in the previous lemma, hence at this point we have bounded uniformly 
all elements in T{di,d 2 , ds, ^4 + A^a^/ 2 - 2 )~^- Observe that in order to derive the bounds for the off-diagonal 
blocks, we used only the bound for the last row of Wu. All elements in the last row of T(di, ^ 2 ; ds, ^4 -|- 
A^ajv/ 2 - 2 )~^ satisfy |r(di, d 2 , da, d 4 -|- A^ajv/ 2 - 2 ) 4 /l ^ A“^ for j = 1,... ,4. It is easy to see that, if we now 
consider r(di, d 2 , ds, d 4 , ds, dg -|- X‘^aN/ 2 - 3 )~^ 1 by a similar argument we obtain for j = 1,..., 6 the bounds 

T(di, d2, da, d^, dg, dg -|- A^ajv/ 2 - a)g/ < A ^ 

T(di, d 2 , da, d^, dg, dg -|- A^ajv/2-3)5/ — 


Finally, from the presentation above follows that assuming that all of the elements in the last row of 
the inverse block T(di,..., d 2 k-i, d 2 k + A^a 7 v/ 2 -fc)~^ are bounded on absolute value by A“^, and the rest by 
22 mA-i, same bounds for the larger inverse block r(di,..., d 2 fc+i, d 2 k +2 + A^ajv/ 2 -(fc+i))~^ will follow. 


thus, we showed that the bounds (77) are propagated for the whole 


-1 


□ 


Lemma 8.4. Let N >0, m> 1, I = 1,... . There exist C{m) > 0 (independent of X and N), sueh that 

for A > 1 the following bounds hold 


E 

E 


{L\ "h 
{L{~\ 


< C{m) ■ {IX) 

< C{m) ■ {IX) 


l+l/2m 


l+l/2m 


for all i = 1,... ,N. 


Proof Here we assume that I is fixed, and we drop the superscript in the notation of T\ aj., 6^, and dj., 
we use simply T, a^, 6^, and dk respectively. Moreover we denote A = IX. 


For the sake of c’ 
From Lemma 


8.2 


z-i 

A • 


-1 


arification let us restrict our attention to the first column of L 

it follows that n = [cA^/^”^] dimensional upper left corner submatrix of ^ is equal to 
T(di,..., dn-i, dn-|-A^a(jv-n)/ 2 )~^) and the absolute values of elements in this matrix are uniformly bounded 
by 2^™'A“^, thus the straightforward estimate for the first part of the sum is 


E 

2=l,...,n 


(ii ) 


il 


< n2^^X~^ = cA^/2”^ • 22™A“^ = 


Next, we are going to show that the terms in remainder X]i=n+i N 
and can be bounded uniformly with respect to the dimension. 

As in the previous lemmas we take the block decomposition of L 


(Li 

i -1 



, i.e. 


obey a geometric decay rate. 


' VL(di,. 

■ ■ ■, dAr)ll 

l^(di,. 

• , dAr)i2 


r(di, . . . , dAr_2) 

-A 

_ W{du. 

■ • , dAr)2l 

l^(di,. 

• , dAr)22 



T{dN-i,dN) 
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From (70) it follows that W{di,... ,dj^) 2 i can be expressed in terms of W{di,... ,djv)ii) namely, for the 
first column of W{di ,..., d 7 v) 2 i (composed of (L(^ )Ar_i,i, {L\ )Ar,i) the identities are 

{L\ )n-i,i = T{dN-i,dN)ii\{L\ )n-2,i-, {L\ )n,i = T{dN-i-,dN)2i'^{L\ )n-2,i- 

Analogously, for {L\ ^) 7 V- 3 ,i, (^a ^)n- 2 ,i we have 

(Aa )n- 3,1 =T{d]Sl-Z,dN- 2 +'^^ai)ll'^{L’‘x )Ar-4,l, (Aa )n- 2,1 = T{dN- 2 ,,dN- 2 +^^ai) 2 i\{L!‘x )iV-4,l- 

From repeating this argument we obtain 

(Aa )N-b,i =T{dN-r)-,dN-i +^ 0 , 2)11 XiLf'x )Ar-6,i5 (Aa )n-a,i = T{dN-b,dN-A +^a2)2i'^{L\ )Af-6,i) 


(-^A )n+l,l — + A"'a(jV-n)/2-l)ll^^(-^A )nl) (-^A )n+ 2 ,l — 7"(dn+l, dn +2 + A^O(7V-n)/2-l)21^) 

This is a recursive series, all elements can be expressed in terms of {V'x ^)ni- Therefore 




I 




E E 


- 1 , 


N-l 


'-'Jh 


i=n+l^...,N 


j=n+l 


where 


Cn+l 

Cn+2 

Cn+3 

Cn+4 


:= T{dn+l-,dn+2 + A^a(jv_„)/2-l)i/A, 

:= T{dn+l-,dn+2 + A^a(jV-n)/2-l)2/A, 

:= T[dn+‘i-,dn+A + A^a(7V-n)/2-2)l/A • r((in+l,dn+2 + Aa(Ar-n)/2-l)21^ A, 
:= T{dn+3^ dn+A + ^‘^0‘(N-n)/2-2)21 ^ ' d^{dn+l,dn+2 + A^a(Af-n)/2-l)2/A, 


CN -2 ■= T{dN-l,d]Xf)iiX ■ T{d]Xf- 3 ,dN -2 + Pai)2i^A- • ■T{dn+l,dn +2 + A^a(jV-n)/2-l)21^A, 
CN-l ■= T{dN-l,dN)21^ • T{dN- 3 ,dN -2 + A^ai)2/A- • ■T{dn+l,dn +2 + A^a(jV-n)/2-l)2/A. 

it follows that > aj > 0, and /(a) = 


From Lemma 


8.1 


X‘^-\-djdj^i+djX‘^a 


is strictly positive for all 


j > 1 (as the derivative is positive, compare (66)). Recall that dj = + (j — 1)^”*, and we have the obvious 

inequality A^ + djdj+i + djX^a < A^ + dj+idj +2 + dj+iA^a for a > 0. Therefore the following inequalities are 
satisfied 


^(dn+j) dn+j+l + A — T{dn+li dn+2 X ®(Af—n)/2— 1)11 

T{dn+jjdn+j+l + X^a(^]\f_(^n+j+l))/2)21 — T{dn+l,dn+2 + X‘^a(N-n)/2-l)2l 


< 


dn+2 + 22-A 


< 


A^ + dn+id,^+2 + 2^™'A(in+i 
A 


A^ T dn+ldn+2 

Now to show the claim about the geometric decay, we take n > 

A . dn+2 + 22™A + {n + lf^ + 2'^^X 


< 


A2 + dn+ldn+2 A^ + dn+ldn+2 + 22"^Adn+l A2 + (^2™ + + („ + l)2m) + 22"^A(/2™ + 

24m+l^ + 22™A 


Similar argument to the one used in (651 shows that 

fra + + l)2m 22«^A 


22™a + 22”^A 
<-^ < 


1 

< —. 


A2 + (P™+ n2™)(l2m + (n +l)2m)+22™A(12 "i + n2”^) a2 + 22"iaA + A2 24™+2A2 + 24"i+iA2 + A2 2A 
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We thus demonstrated that for any N we have 


E 

2=n+l,...,7V 




< 




2E 

2=1 




< 




(78) 


Now taking C{m) > + 2^™’+^ we obtain the claim. 

To conclude, observe that the bound holds for the first row, as the matrices and ^ commute 
(L^ = D + A, where D is a diagonal, A is a skew-symmetric matrix). The bound is true for any other 
column/row, to see this note that for each column there are at most n = [cA^/^”*] elements beyond the 
geometric decay regime, therefore the bound is true for any column of 

□ 


8.1 Proof of Theorem 15.81 

Using lemmas presented in this section we prove the main result with inverse matrix bounds 

Theorem 


5.8 


Let / = 1,... ,iV. Let L^^ be the matrix given by (26), L^ be the matrix {PbL\) ^. 
The following estimates hold for the matriees L^^ ^ (diagonal submatriees of L^ ^). 

\\L{-\\b^,^<2^^{1\)-\ 

||Tr'||p^u<C2(m)(/A)-i+V™. 


The following estimates hold for the matrix Lf 


-1 


Lf <Ci(m)A-'+'/2m^ 

rB~^\\ ^ r)2m\—l 

E ^ A , 

Tr'llzi^d<^2(m)A-i+V-. 


Proof The uniform l^ estimate for each column of Lf ^ follows directly from Lemma 


l°° estimate for each column of Lf ^ follows directly from Lemma 


8.3 


8.4 


The uniform 


In order to estimate uniformly the gradient norm of Lf we are going to consider two cases separately. 
Let c > 2(2™-+^)/^”^, a = LAi/(2m 1) I where [•] is the integer part. Let us demonstrate the result for the 


first column of Lf 


Case I For I < a we split the sum 

N n N 

+ ^ =i: + E + ^ (L-')„ 

^ k=n+l 


k=l 


where n = l^/'^'^X i/[2"^(2m i)]q, — |■g^l/ 2 m_^l/ 2 mj particular choice is due to technical reasons). The 
finite part of the sum above can be estimated 


fc=l 


< max 
k=l,...,N 


n 

^ Y,l + k< 2‘^^{lX)-^ (in + n^) < ci (ZA)”^ ii/m^i/m ^ (;a)-^+^/" 


k=l 
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where we estimated In = ;( 2 m-l)/ 2 m^l/m_^l/ 2 m ^ ^ ^ from 

the proof of Lemma |8.4| it follows that the remaining part of the sum is within the geometric decay regime, 
therefore we can estimate like in (78) 

N 


oo I oo 7 

I: (/+.) (Lr\. < 


k=n+l 


k=l 


k=l 


in the last inequality we used the estimate from Lemma 
The final uniform bound for this case is 

N 


{l + l)<C2\-\ 


E (' + '=) 

k=l 




Case II For I > a. 

For this case we have di > 2^™'+^/A, and therefore from the proof of Lemma 
column is within the geometric decay regime, therefore the whole column can be estimated like in (78) 


8.4 


it follows that the whole 


N 


(z+fc) 


k=l 


V 


-1 


k,l 


< 


(L-) 


1,1 


^ oo 7 r. 


k=l 


k=l 


2k 


< C3 


(L-) 


1,1 


I < C3X 


-1 


Final bound The bound in Case I is clearly of higher order, hence it is the final uniform bound. The 
bound is true for other than the first columns, as there are at most n = elements beyond the 

geometric decay regime. 


□ 
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